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摘 要 :边界 条 件 为 指数 衰减 防 数 的 一 维 热传导 问题 ,Laplace 正 变换 后 泛 定 方程 的 通 解 与 边界 条 
件 象 函数 的 乘积 组 合 求 逆 困 难 ,问题 难以 求解 。 应 用 Laplace 变换 性 质 ,将 边界 条 件 用 算 符 运行 于 
变换 过 程 中 ,建立 一 维 热传导 问题 通用 的 理论 解 ;再 将 特定 问题 的 边界 条 件 代入 通用 理论 解 ,获得 
相应 的 解 。 根 据 边 界 条 件 为 指数 衰减 函数 时 的 解 ,讨论 其 与 经 典 问题 解 的 转换 关系 ,并 结合 实例 ， 
利用 测 点 温度 变化 曲线 的 拐点 求 算 模型 参数 。 文 中 提出 的 求解 方法 ,边界 条 件 不 参与 繁杂 的 交换 
) 过 程 ,解法 便捷 。 
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The solution to one-dimensional heat conduction problem bounded 
by the exponential decay condition and its application 
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Abstract : For the one-dimensional heat conduction problem bounded by the exponential decay condition , it 
is hard to find the inversion of the product of the general solution of the universal equation and the image 
function of the boundary condition after Laplace transform, thus increasing the difficulty in solving the 
problem. The general theoretical solution of one-dimensional heat conduction problem is established by u- 
sing properties of Laplace transform and operating the boundary condition during the transformation process 
with operators. For specific problems, the boundary conditions are substituted into the general theoretical 
solution to obtain the corresponding solutions. According to the solution when the boundary condition is ex- 
ponential decay function ,the transformation relationship between the solution and that of classical problem 
is discussed , and the corresponding parameter for the model is calculated by using the inflection point of 
the temperature measurement variation curve. The proposed solving method is convenient because the 


boundary conditions do not participate in the complicated transformation process. 
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热传导 方程 ,是 偏 微分 方程 中 最 重要 的 发 展 方 
程 ,是 空气 动力 学 、 弹 性 力学 、 电 动力 学 .流体 动力 
学 、 近 场 动力 学 .量子 力学 ,等 离子 体 物理 数字 信号 
处 理 等 领域 的 基本 理论 方程 ;热传导 问题 的 研究 
方法 ,对 上 述 领域 相关 问题 研究 大 多 有 借鉴 意义 ; 
也 即 ,热传导 问题 的 研究 成 果 可 为 上 述 力学 领域 的 
研究 提供 理论 支撑 。 

对 于 一 类 边界 作用 下 的 温度 场 ,将 边界 处 的 温 
度 7(1) 假 定 为 已 知 常数 A7,( 即 ,边界 温度 在 初始 
瞬时 获得 升幅 A7, 并 保持 不 变 ) , 则 构成 经 典 的 一 维 
热传导 问题 ““ ;此 时 ,可 利用 Laplace 变换 等 积分 变 
换 方 法 求解 ,这 也 是 这 类 问题 的 基本 解 。 
AN 自 Joseph Fourier 于 1822 年 建立 起 第 一 个 热 传 

TRUY R ,热传导 问题 一 直 是 数学 物理 方法 研 
究 画 基础 内 容 之 一 ,但 时 至 今日 仍 有 一 些 比较 党 
让 乓 问题 难以 直接 应 用 Laplace , Fourier 等 常用 的 变 
换 对 法 求解 ,如 边界 温度 瞬时 上 升 AT 后 按 牛顿 准 
却 定律 衰减 ”的 传 热 问题 污染 源 浓 度 持续 衰减 下 
BUSES 等 问题 。 当 边界 条 件 o Co) 为 指数 衰减 
HERY AT e “时 ,利用 Laplace 变换 进行 模型 求解 
过 程 中 ,需要 对 模型 正 变换 后 常 微分 方程 的 通 解 
ex 瑟 -xVs/a) 与 边界 函数 的 象 函 数 1/(s + A 的 乘 
FUHR expl- x/s/a) / Gs +A) 进行 求 逆 , 而 这 一 乘 
积 组 合 求 逆 困 难 , 因 此 无 法 直接 求解 。 
斧 由 于 复杂 函数 边界 作用 下 的 热传导 问题 求解 析 
解 过 程 的 复杂 性 ,导致 这 方面 的 理论 研究 和 应 用 较 
2p? ;但 解析 解 在 相应 研究 中 又 显得 十 分 重 
要 中 。 对 上 述 问题 ,直接 应 用 现 有 文献 无 法 求解 ， 
故 本 研究 根据 Laplace 变换 的 特点 ,对 边界 条 件 函 数 
不 直接 进行 变换 运算 ,而 是 用 算 符 将 其 运行 于 求解 
过 程 中 ,系统 应 用 Laplace 逆 变 换 的 卷 积 定理 、 微 分 
性 质 、 卷 积 微分 “7 ,建立 问题 的 通用 理论 解 ,再 将 边 
界 条 件 代入 通用 理论 解 , 求 出 实际 问题 的 解析 解 ; 这 
种 求解 方法 ,可 为 类 似 问题 求解 提供 一 种 不 需要 繁 
杂 变换 过 程 的 便捷 求解 方法 。 


1 基本 模型 


以 一 端 设 置 有 热源 的 薄 层 材料 为 例 (图 1)。 材 
料 中 传 热 问题 的 特点 ,可 概括 为 :试验 材料 均 质 各 向 


同性 ,在 x 方向 上 无 限 延 伸 ; 材 料 一 端 设置 成 有 形 
Dirichlet 边界 的 热源 ,材料 与 边界 外 表面 代为 隔 热 
面 ;材料 与 边界 的 初始 温度 T(x,0) =0 ,温度 场 过 余 
温度 记 为 7T(%,t) ;边界 温度 随时 间 的 变化 函数 记 为 
e) ;热源 在 薄 层 材料 中 的 传 热 为 一 维 热传导 ,边界 
附近 温度 场 变化 见 图 2。 


热源 


d 
图 1 试验 材料 示意 图 
Fig.1 Schematic diagram of the tested material 
7 


T(x,t) 


T(x,0) 


0 x 
图 2 Dirichlet 边界 附近 温度 场 空间 变化 示意 图 
Fig.2 Spatial variation diagram of the temperature field 


near the Dirichlet boundary 


上 述 热 传导 问题 ,可 写成 数学 模型 (I ) 如 下 。 


CR ugue +0, t>0 (1) 

ot Ox 

T(x,t) |,- 2 T(x,0), x >0 (2) 
T(x,t) |,., 2 T(0,0) * g(1), t>0 (3) 

T(x,1) | „n 2 T(x,0), t>0 (4) 


其 中 a 为 固体 材料 的 热 扩 散 率 或 导 温 系数 ,m/s。 
4 u(x,t) = T(x,t) -7T(x,0), 则 模型 ( 工 ) 可 
RAHIRA), B 


ðu Fu 
一 =4 7—5, 0<x< *o,1»0 (5) 

ðt Ox 
u(x,t) |,_o =0,x>0 (6) 
u(x,t) |0 =p), 1>0 (7) 
u(x,t) sue =0,1>0 (8) 


在 g(t) 有 确定 函数 如 经 典 模 型 中 的 常数 AT, 
时 ,模型 ( 卫 ) 可 利用 Laplace 变换 进行 求解 。 


2 通用 理论 解 


在 不 依赖 p( 的 积分 变换 与 逆 变 换 过 程 的 前 
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提 下 ,根据 变换 的 性 质 , 先 给 
论 解 。 

对 模型 ( 工 ) 求 对 于 + 上 的 Laplace 变换 ,可 得 模型 
( 亚 ) 如 下 。 


全 出 这 一 类 模型 通用 的 理 


-0 (9) 
ul. =Lle(t)] (10) 
u|, ,=0 (11) 


其 中 :为 关于 的 Laplace 变换 的 象 函 数 ;s 为 

Laplace 算 子 ;LL 和 ”分别 为 Laplace 变换 算 符 与 逆 
变换 算 符 。 

变换 过 程 中 ,边界 条 件 式 (7) 变换 为 式 (10) 时 ， 

不 直接 求 其 象 函 数 , 而 只 用 算 符 将 其 运行 于 过 程 中 。 

_ 模型 ( 亚 ) 中 式 (9 ) 的 通 解 为 


qe -oep[ > x) E TEE x] 


SN (12) 
其 由 ， cl、\c 为 待定 常数 。 

@ 由 定 解 条 件 式 (10) ~ (11) ,模型 ( 亚 ) 的 特定 
解 为 


ü(x,s) =L[ g(t)] E NES Z (13) 


Okeh, 4 TO) 是 指数 衰减 型 函数 AT en" 
MEELTCO] -L(AT,e") = AT,/(s +A), X F 
eipt- xV/s/a) / (s * A)3K Laplace 道 变换 困难 ; 所 以 ， 
对 C&T,e “不 直接 求 其 象 函 数 的 条 件 下 ,利用 Laplace 
变 医 的 性 质 ,用 算 符 将 边界 函数 运行 于 变换 过 程 中 。 
己 对 式 (13) 进行 Laplace 逆 变 换 , 注 意 “Laplace 道 
ain e pun" ,有 

u(x,t) =L '[u(x,s)] 


-Ze(D)) - ew[ TER x)| 
=L '[L(¢(1))] =L [e| TERR] 
-e L^ [ew HERE (14) 


其 中 * 为 卷 积 算 符 。 
注意 应 用 “互补 误差 函数 erfe(z) "的 Laplace 逆 
变换 , 则 有 


eT net 


2 +0 -2 
m o " e d ( 15) 
ES 
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为 将 式 (15) 应 用 于 式 (14) ,和 需 对 式 (14) 进 一 步 
整理 展开 ,由 Laplace 道 变换 的 “微分 性 质 ”, 有 


a r9 TER jJ 
PES (16) 
由 式 (14) 和 式 (16) ,有 
R EE «| [e ( (17) 


根据 Laplace 变换 的 “ 卷 积 微分 ”性质 可 得 


"MO «ale (2 +ọ(t) - LE 


|= 


t=0 


ees] * Hewl , or) B . erfe( 7. 
(18) 


Herte 了 全 =0, 由 式 (17) ~ (18) 整 理 


xa). 
可 得 


u(x,t) 2g(t) lios . ee zyz] * 
— \ dlo()] 
eer * ? (19) 
注意 u(x,t) 2 T(x,t) - T(x,0) ,应 用 卷 积 的 交 
换 律 , 式 (19) 写 成 积分 形式 为 


Tx d) = TO D) EE |. erte 7 


十 


[ H Wie cm m 5) 


式 (20) 是 在 对 p(t) 未 进行 直接 变换 的 条 件 下 ， 
获得 的 模型 解 , 解 中 含有 e CO) ;具体 到 实际 问题 应 
用 中 时 ,还 需要 代入 已 知 的 g() 后 再 进一步 展开 ， 
方 能 获得 实际 问题 的 解 。 也 即 , 式 (20) 对 任意 9g(?) 
都 成 立 , 故 式 (20) 是 模型 的 通用 理论 解 。 

需要 说 明 的 是 ,上 述 方法 求解 过 程 中 ,虽然 p(1) 
未 进行 形式 上 的 变换 ,但 pg(i) 实 质 上 参与 了 变换 与 
逆 变 换 过 程 ;所 以 ,p(t) 必须 满足 Laplace 变换 对 函 
数 的 基本 要 求 , 即 : 函 数 g(1) 在 有 限 区 间 内 可 积 "; 
在 上 趋 于 无 穷 大 的 时 候 ,1p(2) 1 与 衰减 因子 的 乘积 


d — (20) 
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ET 0; 指数 衰减 函数 AToe “为 边界 条 件 作用 下 的 
传 热 问题 ,显然 满足 上 述 要 求 。 


3 边界 为 指数 衰减 函数 时 的 解 


3.1 问题 的 解 


边界 温度 瞬时 升 高 AT 后 按 和 牛顿 冷却 定律 自然 
冷却 , 即 边界 温度 e(t) 2 T(0,0) +AT, e", HF 
T(x,1) 为 温度 场 过 余 温度 , 记 7(0,0) 2 T(x,0) =0， 
M o) =AToe“, 代 入 到 式 (20) 中 ,有 


T(x,t) - Af, veter) E 
A- AT, [ew ; ete cn) (21) 


CN 将 边界 条 件 直接 代入 这 类 问题 的 通用 理论 解 式 
COEM ES UIT NE TII TE MS 
aEBEEHEPR EE WR R ZOL RANE F DU P IE, AH 
PAS UE HEC JE R PRU A [s] RR fe te T 
参考 。 


32 EM 


IU A S0 时 ,此 条 件 下 p(1) = ATu, 对 应 的 物理 
意义 是 边界 温度 瞬时 变化 ATu 后 保持 不 变 ,这 也 是 
经 旺 问 题 中 的 边界 条 件 。 

所 此 时 , 式 (21) 右 端 第 二 项 为 0, 转化 为 经 典 问题 
的 十 , 即 经 典 问题 的 解 是 式 (21 ) 的 特例 。 


4 ” 解 的 应 用 


4.1 模型 参数 求 算 方法 


在 试验 中 ,利用 温度 场 温 度 动态 监测 数据 ,计算 
试验 材料 的 热 物 性 参数 .也 即 模型 中 的 参数 热 扩散 
率 或 导 温 系数 ,是 这 类 问题 研究 的 重要 目的 之 一 。 

借鉴 文献 [16-18 ] (河渠 附近 潜水 渗流 力学 模 
型 ) 提 出 的 方法 原理 ,距离 边界 x 处 ,温度 变动 速度 
w(x,t) =97T(x,t)/9i, 由 式 (21) 得 


itd uei ve -£j +AAT erte [7 
(22) 


则 yy(xw,?) 的 变化 速率 为 
Op (x,t) _ . 1 eme " x. 
ai = AT, uw e ?| ) 
| -2 +Z a] (23) 
由 式 (22 ) 可知 ,ywy(x,t)-t 曲线 上 ,存在 一 个 
拐点 。 
令 9w(x,t)/9t=0, 求 得 拐点 处 的 时 间 c, 2 


ED 
| a 


根据 距 边界 x 处 的 材料 温度 动态 , 绘 出 (x,t)-t 
曲线 (此 时 x 为 一 确定 值 ) ;利用 曲线 的 拐点 ,以 及 
边界 温度 实际 变动 数据 而 获得 的 AT,、A( 其 中 ,A >0 
表示 边界 温度 下 降 ;A <0 表示 边界 温度 上 升 ) ,计算 
出 模型 参数 a。 

值得 指出 的 是 ,不 论 A 的 正 负 ,在 AT A 为 固 
定 值 的 单调 函数 型 边界 条 件 作 用 下 ,温度 场 中 任意 
一 点 的 温度 变化 速度 曲线 不 可 能 存在 2 个 拐点 ; 即 
关于 六 算式 (24) (25) 只 有 其 中 一 个 是 合理 的 。 妆 
A «0 时 ,算式 (24) 中 i。 <0, 与 问题 的 物理 意义 不 
符 , 即 式 (24) 不 具 普 适 性 ; 式 (25 ) 的 计算 结果 不 产 
生 上 述 矛 盾 , 即 式 (25) 具 有 普 适 性 。 故 可 依据 温度 
场 中 x 点 实测 y(x,i)-i 曲线 ,由 式 (25) 计 算出 模型 
参数 a (IERT, AT, A w 均 已 知 ) 。 


4.2. 实例 研究 


在 地 源 热泵 应 用 中 ,地 下 埋 管 换 热 器 的 设计 万 
为 重要 ,而 其 设计 的 关键 要 素 就 是 获得 准确 的 岩 土 
热 物 性 参数 ”] 。 本 研究 在 合肥 工业 大 学 建筑 艺术 
馆 地 源 热泵 示范 工程 项 目的 5 号 观测 孔 取 一 岩 土 试 
件 , 对 岩 土 的 热 扩 散 率 进行 相应 的 试验 。 该 试 件 
(图 1) d=3.0mb=1.5 m,c 20.3 m; Æ b 的 一 端 
预 设置 一 外 径 为 0.25 m 的 钢管 (钢管 及 岩 土 试 件 四 
AREF) ,进行 隔 热 保护 。 在 试 件 中 距 钢 管 0.3 m 
处 设置 测 温 光 纤 , 对 温度 进行 连续 测试 。 

在 一 次 连续 10 的 试验 中 ,混凝土 试 件 的 初始 
温度 与 室内 空气 温度 都 为 18 ;试验 开始 时 ,将 
36% 的 油 迅 速 注 入 钢管 中 , 之 后 油 温 缓慢 自然 冷 
却 ,至 试验 结束 时 油 温 为 22 CC 。 试 验 中 ,充满 油 的 
钢管 形成 岩 土 试 件 的 边界 , 油 温 自 然 冷 却 符合 牛顿 
冷却 定律 ; 距 边 界 0.3 m(x 20.3 mm) 处 , 测 点 温度 
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T(x,t) 及 温度 变动 速度 (x,t) ,摘录 整理 见 表 1、 
图 3。 

由 图 3,y(x,t)-t 曲线 拐点 出 现时 间 1, =2. 65 h; 
利用 式 (25 ) , 求 得 a =1. 94 x 10 mm/s, 这 在 试验 地 
段 岩 性 热 物 性 参数 的 经 验 值 范围 内 汪 ) 。 

表 1 x=0.3m 处 温度 测量 值 


Tab. 1 Temperature measurements at x =0.3 m 


t/h T(x,t) /*C y(x,t) /(C * h^!) 
0.5 18.14 0.28 
1.0 19.35 2.42 
1.5 21.19 3.68 
.0 23.33 4.28 
B. 25.61 4.56 
.0 27.85 4.48 
= 29.92 4.14 
.0 31.85 3.86 
.0 35.21 3.36 
.0 38.20 2.99 


D(C «a 


v 


O wu) 
—Er- T(x,t) 


图 3 x=0.3 m 处 温度 变化 速度 曲线 


Fig.3 Temperature variation velocity 


curve at x =0.3 m 


5 tà it 
通过 上 述 研究 得 出 结论 如 下 。 

1) 对 于 边界 条 件 为 指数 衰减 函数 AT e" D 
维 热传导 问题 ,往往 难以 直接 利用 Laplace 变换 进行 
模型 求解 ;此 时 ,系统 应 用 Laplace 道 变换 的 卷 积 定 
理 、 微 分 性 质 、 卷 积 微分 性 质 ,建立 模型 的 通用 理论 
解 ,再 将 边界 条 件 代 入 通用 理论 解 , 求 出 实际 问题 的 
解析 解 。 

2) 本 研究 提出 的 通用 理论 解 ,对 一 维 热传导 问 
题 具 有 通用 性 ;这 将 为 复杂 边界 条 件 下 的 类 似 问题 
研究 ,提供 一 种 省 略 繁杂 变换 过 程 的 简单 .便捷 的 求 
解 方法 。 
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